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Extending the Minimal Supersymmetric Standard Model (MSSM) to explain small neutrino 
masses via the inverse seesaw mechanism can lead to a new light supersymmetric scalar partner 
which can play the role of inelastic dark matter (iDM) . It is a linear combination of the superpart- 
ners of the neutral fermions in the theory (the light left-handed neutrino and two heavy Standard 
Model singlet neutrinos) which can be very light with mass in ~ 5 — 20 GeV range, as suggested 
by some current direct detection experiments. The iDM in this class of models has keV-scale mass 
splitting, which is intimately connected to the small Majorana masses of neutrinos. We predict the 
differential scattering rate and annual modulation of the iDM signal which can be testable at future 
Germanium- and Xenon-based detectors. 
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Introduction - A plethora of cosmological observations 
have clearly established the existence of a dark matter 
(DM) component to the matter budget of the Universe, 
which is about five times that of the visible matter con- 
tribution. The mass and interaction properties of the 
DM are however not known, and experimental efforts are 
under way to detect it via its scattering against different 
kinds of nuclei, which will not only provide additional 
direct evidence for its existence, but will also reveal the 
detailed nature of its interaction with matter. Since no 
particle in the Standard Model (SM) can play the role 
of DM, this will also be a sensitive probe of physics be- 
yond SM and supplement the new physics search at Large 
Hadron Collider (LHC). Supersymmetric extensions of 
SM provide one class of such models with several nat- 
ural candidates for DM (e.g. neutralino, gravitino, etc) 
when i?-parity symmetry is assumed. However, if the 
DM mass turns out to be in the few GeV range, as sug- 
gested by some recent experiments the minimal ver- 
sion (MSSM) may need to be extended to accommodate 
this, mainly due to its inability to reproduce the observed 
DM relic density 0. 

Another reason for considering extensions of MSSM 
is to explain small neutrino masses. It would be inter- 
esting to see if the same extensions can also provide a 
new DM candidate and determine its properties. Simple 
ways to understand the smallness of neutrino masses are 
by adding one or more heavy SM singlet fermions to the 
MSSM . In these cases, the superpartner(s) of the sin- 
glet neutrino(s) with a small admixture of the left sneu- 
trino can play the role of DM. In this paper, we focus 
on the supersymmetric inverse seesaw models Q where 
one adds two SM singlet fermions N and S. In these 
models, there are three lepton number carrying electri- 
cally neutral fermions per family, namely (y,N,S). The 
DM particle is the lightest super-partner (LSP) of the 



model which can be a linear combination of the super- 
partners of (y, N, S) 0-0| . Current literature on the sub- 
ject discusses two classes of such models. In the first 
class, inverse seesaw is considered within the framework 
of MSSM however, in these models, one needs to 
omit terms in the superpotential that are allowed by 
the symmetries of the Lagrangian. The second class 
of inverse seesaw DM models extend the gauge symme- 
try to SU(2) L x U(1)y x U(1)b-l so that the inverse 
seesaw mass matrix arises from a B — L gauge symme- 
try However, the B — L gauge symmetry discussed 
in Ref. [f| does not arise from a grand unified theory 
(GUT). Yet, a third class uses global B — L symmetry to 
restrict the inverse seesaw matrix to the desired form [f| . 
In this letter, we extend the MSSM gauge group to 
the supersymmetric Left-Right (SUSYLR) gauge group 
SU(2) L x SU(2) R x U(1) B -l so that not only does the 
inverse seesaw matrix arise naturally but, as was shown 
in Ref. Q, this model can emerge as a TeV scale theory 
from 50(10) GUT implying different dynamical proper- 
ties of the DM particle than the previous woks. 

Working within this new class of models, we find the 
following results: (i) A linear combination of the super- 
partners of the new SM singlet fermions can be the LSP 
with mass from a few GeV up to about 100 GeV without 
running into conflict with known low energy observations, 
(ii) The S'-fermion, which is given a small lepton num- 
ber violating mass to understand small neutrino masses, 
leads to a splitting of the above complex scalar LSP into 
two closely-spaced real scalar fields, the lighter of which 
(we'll denote it by xi) is the true DM field and is ac- 
companied by its slightly heavier partner field (x%) with 
a mass difference of order keVs. A consequence of this is 
that the direct detection process involves a dominantly 
inelastic scattering mode with the nucleus (A/") where 
Xi + N{A, Z) -> X2 +J\f(A, Z) [§], and can therefore be 
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tested in direct DM detection experiments [ioj . The in- 
elastic property arises naturally since the gauge Noether 
current coupling to the Z (and Z' in models with ex- 
tended gauge symmetries) necessarily connects xi to X2] 
also, any possible elastic contribution (mostly through 
the Higgs mediation) is highly suppressed due to small 
Yukawa couplings to light quarks. We believe this is an 
important point which has not been properly emphasized 
in earlier papers, (iii) The new TeV scale gauge dynam- 
ics in these models leads to new annihilation mechanisms 
for DM in the early Universe responsible for its current 
relic density. 

General Structure of Super ■symmetric Inverse seesaw - 
The inverse seesaw model [3( for neutrino masses uses the 
following mass matrix involving the (v, N, S) fields: 
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where we have suppressed the family in- 
dex. This leads to the neutrino mass formula 
M v = MdM^hs ( m dM^ 1 ) T . Here the smallness 
of the neutrino masses arises from the small parameter 
fiSt which for Mr ~ TeV and mo ~ GeV, has a value 
in the keV range. 

Neglecting the keV-scale lepton number violation ef- 
fect, the mass eigenstates are complex scalars in the basis 
of (V, W,S), and the lightest sneutrino eigenstate can be 
written as 



X 



i = + {U f )i Ni Nj + (U*)i Sl Si (2) 



-Msusy as in some region of parameter space in SUSYLR, 
the splitting can be enhanced. 

The mass matrix in Eq. ([1} arises in the SUSYLR 
model Q, after symmetry breaking, from the superpo- 
tential 

W = Wmssm + Y V L<5>L C + Y S S4> R L C + -SfisS , (4) 

where $ is a bi-doublet under SU(2) L x SU(2) R with 
B — L charge zero, and (j)R is a right-handed doublet 
with a B — L = 1 responsible for B — L breaking. 

Relic Density - From Eq. @, we see that the sneu- 
trino DM is a linear combination of the (v, N>, S) fields. 
Therefore, as shown in Fig. [1] its annihilation channels 
involve three major contributions, one from each compo- 
nent. Note that the second and third channels are the 
new contributions in SUSYLR. To leading order, the ex- 
pressions for the annihilation cross sections in the s- and 
i-channel are respectively given by 
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where U is a 9 x 9 unitary matrix that diagonalizes the 
full neutrino mass matrix given by Eq. ([1]) . We note here 
that since the entries in the N, S sector of the sneutrino 
mass matrix are ~ TeV, to get an LSP in the GeV-rangc, 
we do need some fine tuning; however, as shown in the 
next section, we are able to get the right relic abundance 
even with such low-mass DM. Also, from universality ar- 
guments, the model requires the sneutrino LSP to be 
always below ~ 100 GeV, beyond which the lightest neu- 
tralino becomes the LSP. 

When lepton number violation is invoked, the split- 
ting terms J2m,n=i ^mnXmXn can be generated in the 
sneutrino sector (similar to those in Ref. l_l|), and up to 
leading order in the lepton number violating mass term 
Us, the mass splitting of LSP can be written as 



where K f = (% - Q/sin 2 6 W ) 2 + (Q/ sin 2 d w ) 2 \ N c = 3(1) 
for quarks (leptons), v = yfl — 4A/ 2 / s is the speed of 
the DM particle in the center-of-mass frame, and C(o,i,2) = 
y^ 3 _ 1 \U( Vt N,s)n \ 2 ■ W e note that both s- and t-channel anni- 
hilations in our case are p-wave scattering, as expected from 
symmetry arguments. For low- mass DM (M x < 20 GeV), 
and assuming Ys ~ 0(1) and M~ ^ 500 GeV, the f-channel 
involving only leptonic final states turns out to be the domi- 
nant contribution in our case. For this reason, we do not show 
the interference term between s- and t-channels for leptonic 
final states in Eq. ([5]). 





SM X = A\A 11 \/M X 



(3) 



FIG. 1. The dominant annihilation channels of the sneutrino 
DM in SUSYLR model. 



where generically An ~ /i^MsuSY- If M x is also of order 
of the SUSY breaking scale (assumed to be around TeV), 
the splitting is of order /is, and if M x is much lower than 



The annihilation cross section for the Z'-channel is 
suppressed compared to the Z-channel by a factor 
{ci/c ) 2 {M z /M z ,) 4 . Also, we find that the correct DM relic 
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FIG. 2. In the left-panel, the purple (black) points correspond 
to the allowed values in the Com -relic density plane, for the 
mixed sneutrino DM in our model; the vertical line shows 
the upper limit for Co from invisible Z-width constraint. The 
right-panel shows the scatter plot of relic density prediction 
for light LSP mass. The horizontal shaded region is the 3<r 
limit obtained from 7-year WMAP data. 



tan#jv = Mn/Mp, Mn lP being the DM scattering ampli- 
tudes off neutron and proton respectively f (|q|) is the 
nuclear form factor and /(v) is the velocity distribution of the 
local galaxy. « m in is the minimal velocity needed to generate 
the nuclear recoil energy E r , which can be written as [f| 
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where <5 is the mass gap, and 8 = corresponds to the case of 
elastic scattering. 

In the case of elastic scattering conducted by Higgs, gm 
in Eq. ([7]) is the total scattering cross section which can be 
written as 
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density is obtained only for Co < 0.16, as shown in the left- 
panel of Fig. [5] This also agrees with the invisible Z-decay 
width constraint (as shown by the vertical line in Fig. 
The right panel of Fig[2] shows the the scatter plot of the pre- 
dictions for the LSP relic density; we find enough parameter 
range where the correct relic density is reproduced for a light 
DM. 

Direct Detection - As the sneutrino DM (xi) in inverse 
seesaw is a real scalar field accompanied by its slightly heavier 
partner field, it has both elastic and inelastic interaction with 
nuclei in direct detection experiments. 

The direct detection channel mediated by the SM-like Higgs 
boson is due to the interaction term Aft-oXiXij where A is 
mainly the D-term contribution which can be simplified to 
(<?Il c o + ffI« c i) l 'wk/4 for SUSYLR assuming the decoupling 
limit of the MSSM Higgs mixing angles, i.e. large tan/3 and 
a ~ f3 — 7r/2 [12J. After invoking the lepton number violation 
effect, a mass splitting is generated between xi an( A X2, an d 
the interaction term can be rewritten as -|/io(Xi + xi) which 
is clearly an elastic interaction. 

The direct detection channel conducted by gauge bosons 
can be written as 
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where 8w is the Weinberg angle. After the lepton number 
violation term is included, the interaction term is of the form 
iZ M (xid M X2 — X2<9/jXi)- Therefore, the collisions between xi 
and nucleus conducted by gauge bosons is inelastic. 

The differential scattering rate of DM particle on target 
nucleus in direct detection experiment can be written as 
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where E r is the nuclear recoil energy, p xi is the local mass 
density of DM, M x is the mass of the DM particle, on is 
the DM-nucleon cross section, and fi X N is the reduced mass 
of DM and the target nucleus, a and A e $ are defined as 
a = (<t p + o- n )/2 and A cB = S ieisotopes 2ri[Zcos6 N + (At - 
Z) sin#jv] 2 , where Vi are relative abundances of isotopes, and 



whereas in the case of inelastic scattering, crjv can be written 
as 
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where the first term in the bracket is induced by Z boson 
whereas the second term is due to the Z' boson in SUSYLR 
and the last term is due to the interference between the two. 
The factors k p — (f — sin 2 6W) and k„ = (|) are due 
to the different coupling of the vector bosons to proton and 
neutron respectively. 

It is important to note here that for large fractions of the 
left sneutrino component (cq ii 10 -3 ), the scattering is mostly 
dominated by the Z-exchange, and is hence inelastic. The Z' 
contribution to the inelastic channel is always suppressed by 
its mass, and similarly, the elastic cross section is 4-5 orders of 
magnitude smaller than the Z-dominated inelastic contribu- 
tion because the coupling of the Higgs to nucleon is suppressed 
by light quark masses. However, for Co ^ 10 -4 and large ci, 
the inelastic contribution, suppressed by the Z'-mass, could 
become comparable to the elastic counterpart. This is shown 
in Fig. [3] where we have plotted the cross section as a func- 
tion of the DM mass for various values of Co and c\. One 
can see that for small Co and large c\ values (blue and orange 
curves), the cross section is dominated by the elastic channel 
(thin lines) for low mass DM and by inelastic channel medi- 
ated by Z' for M x <; 10 GeV (thick lines), whereas for small 
or zero ci component (pink and green curves), the scatter- 
ing is always dominated by the inelastic channel mediated by 
Z. It is interesting to note here that the current XENON100 
data constrains most of the parameter space of the model 
and for large Co, puts an upper bound on the DM mass. This 
can be seen clearly from Fig. [3] where we have shown the 
XENON 100 limits on scattering cross section for various val- 
ues of the mass gap, starting from zero on the left (red solid 
curve, corresponding to the elastic case) to 8 = 30, 60, 90 
and 120 keV cases. We note that for small mass gaps of order 
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a few keV (as expected in this model), the XENON 100 con- 
straints leave only the low mass iDM (below 20 GeV) open 
in this model. Similar constraints on the iDM mass can be 
obtained using the existing data from other direct detection 
experiments [lfj], but the XENON 100 constraints are found 
to be the most stringent. 
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energy region. Also from the solid blue curves in the first and 
third plots, one can see that the energy regions for large re- 
coil energy and large modulation can be separated from each 
other, which provides a chance to fit the anomaly observed by 
the CoGeNT experiment. In these plots, Aq and A\ are the 
zeroth and first Fourier modes of the differential scattering 
rate; for a detailed definition, see Ref. [1 81 ] . 
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FIG. 3. The model prediction for the scattering cross section 
of the sneutrino DM off a nucleon as a function of its mass, 
for various choices of the left and right sneutrino components. 
Also shown are the XENON 100 constraints for mass gap S = 
(0, 30, 60, 90, 120) keV which were obtained by Feldman- 
Cousins method [l4[ using the 100 days of live data [l5[ . 

The nuclear recoil energy in both elastic and inelastic scat- 
tering has a maximum determined by the escape velocity of 
DM in the local galaxy, and a minimum determined by the 
mass gap for inelastic case. Therefore, the topology of the dif- 
ferential scattering rate for inelastic scattering is very different 
from the elastic scattering, which can be used to determine 
whether the DM is inelastic or not. Furthermore, for inelastic 
scattering, if the mass gap is comparable to the kinetic energy 
of DM, for certain nuclear recoil energy, due to Eq. ((HJ, the 
required velocity is pushed to the tail end of the velocity dis- 
tribution where the motion of the earth has a larger effect and 
therefore the DM annual modulation signal gets enhanced. 

The predicted normalized differential scattering rate and 
annual modulation for Germanium- and Xenon-based detec- 
tors are shown in Fig. [4] with different choices for param- 
eters of the SUSYLR model. The mass of Z' is taken to 
be 1 TeV in both cases. The red and blue curves are for 
(c ,ci) = (0.001,0.1) and (co, a) = (0.1,0.001), correspond- 
ing to the Z' and Z dominance, respectively, in the inelastic 
scattering between DM and target nucleus. The latter case is 
similar to the MSSM version of inverse seesaw. To translate 
the differential rates to experimental quantities, namely the 
electron equivalent recoil energies in germanium detector and 
the Si signal in xenon detector, we have used the quench- 
ing factor and scintillation efficiencies from CoGeNT [TtJ and 
XENON 100 [ill experiments respectively. One can see from 
the first two plots in Fig. [4] that a peak shows up if the scat- 
tering is dominated by inelastic interactions. Furthermore, 
one can see from the third and fourth plots that for inelastic 
case, the annual modulation can be larger than 100% in some 



FIG. 4. The model prediction for differential scattering rates 
(normalized to one) and annual modulations for germanium 
and xenon detectors. The upper two plots show the normal- 
ized scattering rate for germanium (left) and xenon (right) 
detectors. The lower two plots show the annual modulation 
in the same detectors. Red and blue curves are for (co, ci) = 
(0.001,0.1) and (c ,ci) = (0.1,0.001), respectively, while the 
solid and dashed curves are for (M x ,5) — (lOGeV, 20keV) 
and (50GeV, 60keV), respectively. 

We note here that it is also possible to explain the recent 
CRESST-II data for an iDM with mass splitting 0(100) keV, 
consistent with the null results from XENON100 19]. How- 
ever, it is difficult to do so for DAM A/LIBRA, even though 
there still exists a narrow iDM parameter space that does fit 
DAMA/LIBRA data with CRESST-II ;2Q]. For this reason, 
we did not include the DAMA/LIBRA results in our discus- 
sion on annual modulation. 

Comments - As noted earlier, there exists a domain in the 
parameter space for which the dominant annihilation channel 
of the DM is to lepton-anti-lepton final states. The relative 
branching fractions depend on the masses of the RH neutrinos 
and are somewhat model-dependent. If we take the model in 
Ref. [8J as a guide, we expect them to be of similar order. 
These dominant ly-leptonic annihilation modes can, in princi- 
ple, be important in understanding signals from the galactic 
center, radio filaments as well as WMAP haze [2lj]; however, 
we expect these effects to be suppressed in our model due to 
the p-wave nature of the dominant DM annihilation channels. 

Furthermore, we observe that there exist upper limits on 
the direct DM detection cross section from monojet plus miss- 
ing energy searches in colliders [2^| . However, for light spin 
zero dark matter only weakly interacting with the nucleons, 
the collider limits on the spin-independent cross sections are 
weaker than the direct search bounds. 

Also, depending on the sparticle spectrum prediction of 
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the model, the sneutrino LSP can have distinctive collider 
signatures. In particular, if the gluino is lighter (heavier) 
than the squarks, the sneutrino LSP can have interesting four 
(two) jet+like (opposite)-sign dilepton signals with missing 
energy 23]. A detailed collider simulation of these events for 
our model will be given elsewhere [24fl . 

Conclusion - To summarize, we have shown that the su- 
persymmetric inverse seesaw model for neutrinos naturally 
leads to an inelastic scalar DM. The differential scattering 
rate and annual modulation predicted in these models might 
be tested in future direct detection experiments. The DM par- 
ticle mass is found to be strongly constrained by the current 
XENON bounds, and for keV-scale mass splitting for the real 
scalar LSP states (as required by neutrino mass constraints), 
we find an upper limit of around 20 GeV on the DM mass. 
This is consistent with the model prediction for the sneutrino 
LSP mass which is required to be below ~ 100 GeV from 
universality arguments. 

Therefore, we might be able to identify SUSY inverse see- 
saw if from the collider search, we can confirm that the sneu- 
trino is a long-lived particle, and from direct detection exper- 
iments, we observe an inelastic WIMP from the differential 
scattering rate and the annual modulation. 
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